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Abstract 

The mean mass densities of cosmic dark matter is larger than that 
of baryonic matter by a factor of about 5 in the ACDM universe. 
Therefore, the gravity on large scales should be dominant by the dis- 
tribution of dark matter in the universe. However, a series of ob- 
servations incontrovertibly show that the velocity and density fields 
of baryonic matter are decoupling from underlying dark matter field. 
This paper shows our attemps to unveil the physics behind this puzzle. 
In linear approximation, the dynamics of the baryon fluid is completely 
governed by the gravity of the dark matter. Consequently, the mass 
density field of baryon matter pj,(r,t) will be proportional to that of 
dark matter pdm( r ^)) even though they are different from each other 
initially. In weak and moderate nonlinear regime, the dynamics of the 
baryon fluid can be sketched by Burgers equation. A basic feature of 
the Burgers dynamics is to yield shocks. When the Reynolds number 
is large, the Burgers fluid will be in the state of Burgers turbulence, 
which consists of shocks and complex structures. On the other hand, 
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the collisionless dark matter may not show such shock, but a mul- 
tivalued velocity field. Therefore, the weak and moderate nonlinear 
evolution leads to the IGM-dark matter deviation. Yet, the velocity 
field of Burgers fluid is still irrotational, as gravity is curl-free. In 
fully nonlinear regime, the vorticity of velocity field developed, and 
the cosmic baryonic fluid will no longer be potential, as the dynamics 
of vorticity is independent of gravity and can be self maintained by 
the nonlinearity of hydrodynamics. In this case, the cosmic baryon 
fluid is in the state of fully developed turbulence, which is statistically 
and dynamically decoupling from dark matter. This scenario provides 
a mechanism of cohenent explanation of observations. 

1 Introduction 

The mass field of the universe is dominated by dark matter, and only a tiny 
fraction of cosmic matter is in the form of baryonic particles. Gravitational 
clustering of dark matter is the major imputs of cosmic structure formation. 
Collapsed massive halos of dark matter are believed to be the hosts of various 
observed objects. On the other hand, light-emitting and absorbing objects 
are made up of baryonic matter. 

Most of the baryonic matter in the universe actually is made of inter- 
galactic gas, or intergalactic medium (IGM), mainly consisting of neutral 
and ionized hydrogen, neutral and ionized helium [U El E]- For a long time, 
the mass field of the IGM, Pfi(r), is believed to be proportional to dark mat- 
ter field pdm(r)- The similarity relation p&(r) oc pdm(r) should be hold, at 
least, on scales larger than the thermal diffusion or Jeans length of the IGM. 
In other words, the IGM mass field is considered as a point-by-point tracer 
of the dark matter distribution. Obviously, in this case, the baryon fraction 
p&(r)/[pdm( r ) + p&(r)] is independent of r, and equals to the cosmic mean 
fl b /(fl b + fidm) = 0.17 ± 0.01 [HE], where fl b and f2d m are, respectively, the 
cosmological parameters of baryonic and dark matter 

However, the similarity p&(r) oc pdm( r )j is found to be seriously incon- 
sistent with observation in collapsed objects. The baryon fractions in grav- 
itational bound objects generally are much lower than the cosmic mean. It 
raises the so-called missing baryon problem. The baryon fraction is claimed 
to decrease monotonically with the mass of collapsed halos. For X-ray clus- 
ters and groups, the baryon fractions are smaller than the cosmic mean by a 
factor of 2 - 4 [HI [71 [8] . The baryon fraction of dwarf galaxies can be as low 
as about 1% of the cosmic mean [5]. 

Many works have been carried out to explain the deficit of baryon content 
in gravitational collapsed objects, including feedback due to supernovae [10] 
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or Active Galactic Nuclei(AGN) [IIJ[T2] , which blows out hot gas from halos, 
and pre-heating [131 E], which heats baryon gas in the IGM and slows down 
the accrection of the IGM into halos. However, cosmologically hydrodynamic 
simulations show that these mechanisms fail to produce the baryon fraction 
in halos reveiled by observation [T5J [16l |T7] 

Actually, the baryon-dark matter decoupling probably is due mainly to 
the nonlinear nature of the hydrodynamics of the velocity and mass fields of 
IGM. That is, even without extra heating and cooling, the hydrodynamical 
evolution can already provide mechanisms of separating the baryonic matter 
from dark matter. A possible dynamical mechanism driving the baryonic 
matter deviate from dark matter has been addressed in the pioneer work 
of Shandarin & Zeldovich [18J. As the dark matter particles are assumed 
to be collisionless, their velocities would be multi-valued at the intersection 
of trajectories. However, the baryon matter has a single-valued velocity 
field as a compressible fluid. Discontinuities, such as shocks or complex 
structures, will develop inevitably in the baryon fluid at these intersectiones. 
Shocks would then cause the under-dense of baryonic matter in the post- 
shock area. On the other hand, dark matter is not affected by the shocks. 
Thus, discontiunities might prevent baryon matter to follow dark matter. 
Consequently, the baryon fraction in the collapsed halos would be less than 
cosmic mean |19j . 

Therefore, the nonlinear evolution of IGM would be a major reason of 
the statistical and dynamical decoupling between the cosmic baryon fluid 
and dark matter. We report the results of this mechanism. This paper 
will be arranged as follows. §2 is on the linear theory of evolution of IGM 
perturbations. §3 presents the decoupling in moderate nonlinear regime. §4 
studies the fully nonlinear development of the IGM. It shows the cosmic 
baryon fluid to be turbulent on small scales. §5 gives observation evidences. 
The conclusions will be in §6. 

2 Linear evolution of the perturbations of IGM 
mass field 

2.1 Basic equations 

Let us consider a flat universe having cosmic factor a(t) oc t 2//3 , and domi- 
nated by dark matter. We describe the IGM by a mass density field p(x, t) 
and a peculiar velocity field v(x, t), where x is the comoving coordinate. In 
hydrodynamical description, the equations of the IGM consist of the conti- 
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nuity and the momentum (or the Eular) equations as [2U] 



% + -V-(l + 5)v = (1) 

at a 

^ + ( v . V )v = --Vp-V0 (2) 

where the density perturbation of the IGM 5(x, t) = [px, t) — p]/p, and p is 
the mean density of the IGM. p is the pressure of the IGM. 

The IGM can be considered as a "passive substance" with respect to the 
underlying dark matter. Namely, the evolution of the IGM mass field p(x, t) 
is governed by the gravity of total matter (baryonic plus dark matter). The 
peculiar gravitational potential <p in eq.(2) satisfies 

V 2 = AnGa 2 p t 5 t . (3) 

where the operator V is acting on the comoving coordinate x, and <5i(x, t) = 
[pi(x, t) — pt]/pt is the perturbation of total density field, with mean density 
Pt = 1/QirGt 2 oc a~ 3 . The potential <ft will be zero (or constant) when density 
perturbation 5 t = 0. 

To sketch the gravitational clustering of the IGM, we will not consider 
the details of heating and cooling. Thermal processes generally are local, and 
therefore, it is reasonable to describe the thermal processes by a polytropic 
relation p(x, t) oc p 7 (x, t). Thus eq.(2) can be rewritten as 

dav lk B T V8 _ , , A . 

— - + (v V v = —L-?—- r--V0 4 

dt pm p {1 + 5) v ' 

where the parameter /i is the mean molecular weight of the IGM particles, 
and m p is the proton mass. In this case, we don't need energy equation, and 
the IGM temperature evolves as T oc p 7-1 , or T = T (l + <5) 7_1 . Exactly 
speaking, there is also a noise in the hydrodynamical equations eq.(2) and 
(4), which are related to the viscous of the IGM fluid due to the fluctuation- 
dissipation theorem. Since the IGM is a "passive substance" , this noise would 
be ignorable in searching for the deviation of the IGM from the underlying 
dark matter. 



2.2 Linear solution of growth mode 

The unperturbed solutions of the density and velocity fields of baryonic mat- 
ter are p = [Q b /(Q b + fi dm )]p t oc a~ 3 , and v = 0. Therefore, the linearizion 
of eqs. (1) and (4) yields 

| + iv.v = (5) 

ot a 
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^ = -7tEL VS _ v , (6) 
at firrip 

where the mean temperature T oc p 7-1 oc a -3 ^ 7 " 1 -*. In the Fourier space, we 
have 

^ + ^CM) = -^(M) + ^(M), (8) 

where v(k, £) = zku(k, t), and the Jeans wavenumber kj is given by k 2 = 
(a 2 /t 2 )( y fim p /'ykbT) oc a 37-4 . When 7 = 4/3, fcj is time- independent. Here, 
we ignore the stochastic nature of eqs.(7) and (8), and the randomness of all 
the variables are given by the random initial conditions. 

We consider only the growth mode of the perturbation of dark matter, 
i.e. 5(k,t) oc a. In the case of 7 = 4/3, the solution of eqs.(7) and (8) is [2T] 

* (k - () = 1 + cir<I+ " /6 + C2r< "" /6 (9) 

where e = (1 — Ak 2 /9kjY^ 2 , and constants c\ and C2 depend on the initial 
condition <5(k, 0) and u(k, 0). Therefore, regardless of the initial condition of 
the IGM, after a long evolution we always have 

5(k,t) = 5 t (k,t), iik^:kj. (10) 

This means that the evolution of the IGM will completely follow dark matter, 
and the initial difference between <5(k, t) and S t (k, t) will be forgotten. Eq.(10) 
yields 

Pb(k) = P d m(k), if k < kj. (11) 

The power spectrum of IGM field, P&(fc) has to be the same as that of dark 
matter P^ik). 

The linear solutions of the IGM have also been studied by using different 
assumptions of the IGM thermal processes other than that used in eq.(4) 
1221 1231. A common feature of these solutions is 



5(k, t) = (1 + decaying terms)(^(k, t) + decaying terms if k kj. (12) 

Only the decaying terms are affected by the initial conditions, and therefore, 
the solutions of eqs.(10) and (11) hold in general, regardless of the specific 
assumptions of the IGM thermal processes. 

The physical explanation of the solutions (10) and (12) is straightforward. 
The IGM is only a tiny component in the cosmic mass field and its evolution 
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is completely governed by the gravity of dark matter in the linear regime. 
After the decaying of the initial conditions, the IGM should follows the same 
trajectory as dark matter on scales larger than the Jeans length. 

In linear approximation, we also have the linear relation as follows: 

5(x,t) = -/V-v(x,t), (13) 

where the parameter / depends on dark matter models [24] 



3 Moderate nonlinear regime of the IGM evo- 
lution 

3.1 Burgers' equation for cosmic baryonic fluid 

In the moderate regime, we may use linear relation eq.(13) for S in eq.(4). 
We have then 

— - + v • V v = — V 2 v - V0 14 

at a 

where the coefficient v is given by 

V ^m p {dhiD{t)/dty 1 } 

in which D(t) describes the linear growth behavior. The term with v in 
eq.(14) describes a diffusion which is characterized by the Jeans length kj. 

To sketch the evolution of gravitational clustering, only the growth modes 
are of interest. For the growth modes, the velocity field is irrotational, and 
we can define a velocity potential by 

v = --Vy?. (16) 
a 

Substituting eq.(16) into eq.(14), we have 

|__L (v , f _ii vW . (17) 

Equation (17) is the stochastic-force-driven Burgers' equations or the KPZ 
equations [25] . It is the simplest nonlinear Langevin equation for growth 
modes [25J I2H]- This equation can also be found with a coarse-graining 
approximation [29J. 

Thus, in moderate nonlinear regime, the cosmic baryonic fluid is similar 
to Burgers' fluid. The second term on the l.h.s. results from convection, 
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and the third on the l.h.s. describes relaxation of the clustering by diffusion. 
The term on the r.h.s represents the gravitational potential of dark matter, 
which provides the initially random perturbations. Clustered structures will 
develop from the initial seeds via the competition of the convection and 
diffusion. 

3.2 IGM as a Burgers' fluid 

The Burgers' equation contains two scales: the dissipation length or the 
Jeans length 1/kj and the correlation length r c of the random force, (p. The 
intensity of the random force <fr can be characterized by the density contrast 
of dark matter 5dm- Many theoretical works on Burgers fluids have been 
done [301 EE E21 E31 EI]- They show that turbulence will develop in Burger 
fluid when the following Reynolds number is large [331 E2] 

{kjr c f/\5l m )^ > 1. (18) 

This condition corresponds to the system in which the dark matter field has 
entered the non-linear regime. In the linear regime, 8d m <C 1, and therefore, 
perturbations on scale r c > 1/kj will not cause Burgers' turbulence in the 
IGM. In the moderate nonlinear regime, 8 dm > 1 ; an d Burgers' turbulence 
on scales larger the Jeans length can develop. Since the dark matter field 
becomes nonlinear first on small scales and then on large scales, one can 
expect that Burgers' turbulence in the IGM develops from small to large 
scales. 

When Burgers turbulence develops in the IGM, its velocity and mass 
density fields will be intermittent and contains discontinuities, or shocks. 
The probability distribution function (PDF) of velocity field is long tailed 
[33] . Since dark matter is not affected by Burgers turbulence, the IGM 
velocity field will dynamically decouple from the dark matter field on scales 
larger than the Jeans length once Burgers turbulence is developed. The 
intermittence of the IGM mass field, as a result of the Burgers equation, has 
also been emphasized in [23] . 

Jeans diffusion will also lead to the decoupling between the velocity fields 
of the IGM v and dark matter Vd m - However, the decoupling given by the 
Jeans diffusion is very different from that results from the Burgers turbulence 
shocks. For the former, the distribution v and Av = v(x + r) — v(x — r) is 
symmetric with respect to the transformation Vd m — > — Vd m , i-e., the velocity 
PDF with Av • Vdm/|vdm| > (acceleration in the direction of Vd m /|vd m |) is 
the same as the PDF with Av • Vd m /|vd m | < (deceleration in the direction 
°f v dm/| v dm|)- For Burgers turbulence, the shocks consist of an acceleration 
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ramp followed by a rapid deceleration. Therefore, the IGM is not symmetric 
between the sections accelerating and decelerating. Since the acceleration is 
due to the gravity of dark matter, the acceleration is generally in the direction 
of Vdm- Thus, the IGM velocity field v is asymmetric with respect to the 
transformation Vd m — > — Vd m - The PDF of Av with Av • Vd m /|vd m | > will 
not be the same as the PDF of Av with Av • Vd m /|vd m | < 0. 

If v — 0, Burgers equation (17) is scale free when the random force <fi is 
scale-free. In this case, the velocity field of the Burgers fluid is self-similar 
and the PDF of the velocity difference Av is scale-invariant. 

In summary, if the IGM can be described as a Burgers' fluid, its velocity 
field should show the following properties: 1.) the peculiar velocity of the 
IGM at a given point will generally be lower than that of dark matter at the 
same point, 2.) the PDF of Av r ■ Vd m /|vd m | > will not be the same as the 
PDF of Av. r • Vdm/|vdm| < 0, and 3.) The PDF of Av r is scaling. These 
predictions have been proven by cosmologically hydrodynamic simulation 
samples [36J. 



4 Fully developed nonlinear regime 
4.1 Vorticity 

Burgers turbulence is qualitatively different from turbulence described by the 
Navier-Stokes equations. The latter generally consists of vortices on various 
scales, while the former is a collection of shocks. These features arise because 
for growth modes of the perturbations of velocity fields is potential, and the 
velocity field is irrotational. Moreover, from the condition eq.(18), we see 
that the Burgers turbulence can only develop in dense areas of dark matter 
5 dm > 1. 

However, in fully developed nonlinear regime, the cosmic baryon fluid will 
on longer be potential, and the turbulence of Navier-Stokes fluid can also be 
developed in high as well as low density regions. This point can be seen with 
the vorticity, which is defined by uji = (l/2)(diVj — djVi) or oj = V x v, where 
% = 1, 2, 3. The dynamical equation of the vorticity u of cosmic baryon fluid 
can be derived from Navier-Stokes equations as follows [37] 

= d t tu + -v • VcJ = -(S ■ to — dcu + -4-Vp x Vp — au), (19) 
Dt a a p 2 

where d = diVi is the divergence of the velocity field. Tensor S is called strain 
rate defined as = (l/2)(9jUj + djVi), and the vector [S • co]i = SijUij. 

A remarkable property of eq.(19) is the absence of the term of gravity 
of mass fields, indicating that the gravitational field of dark matter is not 
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the sources of the vorticity. This point is obvious, as gravity is curl-free in 
nature. In this sense, the vorticity is " non-gravitational" . 

In the linear regime of perturbations, both the terms S • u — du and 
Vp x Vp would be zero. Vorticity Q will decay as a -1 , given by the term 
au. Therefore, as expected, in the linear regime the velocity field of baryon 
matter is irrotational. In other words, the linear evolution of the velocity 
field of baryon fluid is fully governed by the gravity. 

However, in the fully developed nonlinear regime, the vorticity will no 
longer be zero. First, it has been shown that at the moderate nonlinear 
regime, Burgers' fluid will inevitably lead to the development of shocks and 
complex structures, which can yield non-zero term Vp x Vp of eq.(19), being 
called baroclinity. Thus, the vorticity is excited, regardless of the gravity 
of dark matter. Physically, once multi-streaming has developed, complex 
structures, like curved shocks, will lead to deviation of the direction of Vp 
from that of Vp. Consequently, the density pressure relation cannot be 
simply given by an single-variable function equation as p = p(p) |38j and the 
baroclinity term will no longer be zero. 

Once vorticity is initiated, eq.(19) shows that the vorticity of baryon ve- 
locity field can be maintained by itself, i.e. the system actually is nonlinearly 
self-excited and self-maintained. The divergence d = djVj is generally neg- 
ative in regions of clustering. Therefore, the term — du of eq.(19) will lead 
to an amplification of vorticity in overdense regions. On the other hand, the 
vorticity will be stretched by the nonlinear term S ■ Q. Therefore, vorticity 
can also develop in low mass density region. This picture is very different 
from the Burgers' turbulence. 

Vortices are considered as a fundamental ingredient of turbulence. The 
fluctuations of the vorticity field is an important measurement of the tur- 
bulence of fluid SO]. Therefore, the dynamics of vorticity provides a 
non-gravitational mechanism of the decoupling between baryon matter from 
dark matter. 

The turbulence of cosmic baryon fluid can also be seen with the spectrum 
tensors $y(k) and f^-(k), which are, respectively, the Fourier counterparts 
of the two-point correlation tensors of velocity (fj(x + r)fj(x)) and vorticity 
(wi( x + r ) w i( x ))- That is 

<Mk) = JL j ( %(x + r )^(x))e-^r (20) 

%(k) = (235 / M x + r)u^))e-^dr. (21) 

respectively, where (...) denotes average over spatial coordinates x. For a 
homogeneous turbulence, we have [39] 
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a,(k) = [Sijk 2 - hk^uik) - fc 2 %(k), (22) 

and hence, 

n«(k) = k 2 ^ u (k). (23) 

The power spectra of velocity and vorticity fields are defined respectively as 

P v (k) = J±Q u Qz)6{\k\-k)dk; P w {k) = Jhl u {k)5{\k\-k)dk. (24) 

Combining eqs. (23) and (24) yields 

P u {k) = k 2 P v (k). (25) 



This relation can be used to measure the developed level of turbulence. If 
the velocity and vorticity fields of a fluid satisfy the relation given by eq.(25), 
it should be in the state of fully developed homogeneous turbulence. Other- 
wise, it would be less developed. The samples of cosmological hydrodynamic 
simulation have showed that the relation eq.(25) is perfectly hold on scales 
less than 3/i _1 Mpc at z ~ [37]. Actually, the dynamics of eq.(19) reveals 
that the cosmic baryonic fluid underwent a scale-free evolution, which leads 
to a state of fully developed turbulence. 



4.2 She-Leveque intermittence 

The velocity field of a fully developed turbulence is highly intermittent, which 
is measured by the structure functions S p (r) defined as 

S p {r) = (6v>) ~ r c *\ (26) 

where 5v r = [(v(x + r) — v(x)) • r/r] and ( p is called intermittent expo- 
nent. Based on dimensional argument of hierarchical evolution, Kolmogorov 
in 1941 [U] predicted that for fully developed turbulence, the intermittent 
exponent is ( p = p/3 on scales of inertial range. Experimental and numer- 
ical results do not, however, support the p/3 law. It should be attributed 
to intermittency, i.e. turbulence field in inertial range is characterized by 
stronger non-Gaussianity on smaller scales. A remarkable development was 
made by She and Leveque (SL hereafter) [12] . They proposed that the 
non-Gaussian behavior of fully developed turbulence is determined by the 
hierarchical structure originated from the Navier-Stokes equation, and the 
p/3 law should be replaced by 

Cp/Cs = [1 - C(l - ^)]p/3 + C(l - /F), (27) 
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which contains only two dimensionless parameter C and (5. C is the Hausdorff 
dimension of the most dissipative structures, and parameter is given by 
the hierarchical evolution. SL formula eq. (27) is in excellent agreement with 
various experiments of turbulence, including turbulence in compressible fluid. 
The SL scaling law of structure function is considered to be universal for 
characterizing the fully developed turbulence. 

With the hydrodynamic simulation sample of the concordance ACDM 
universe, the intermittent exponent of the velocity field of cosmic baryon 
fluid at redshift z = in the scale range from the Jeans length to about 
16 h -1 Mpc is found to be extremely well described by the She-Leveque's 
universal scaling formula eq.(27) with C — 1 and (3 3 = 1/3 [13]. It indicates 
that the dissipative structures are dominated by sheets. These results imply 
once again that the evolution of highly evolved cosmic baryon fluid is similar 
to a fully developed turbulence. 

When a fluid is turbulent, the kinetic energy passes from large to the 
smallest eddies, and finally dissipates into thermal motion. For cosmic bary- 
onic gas, the evolution of vorticity is also hierarchical. The vorticity evolves 
from large scales to small scales, and finally falls into massive halos to form 
structures, including light-emitting objects. This result strongly indicates 
that the highly nonlinear evolution would lead to the cosmic baryon fluid 
reaching a statistically quasi-equilibrium state satisfying the universal scal- 
ing as that of fully developed turbulence. 

In view of this picture, we can say that in the highly nonlinear regime, 
the statistical properties, especially the intermittent behavior, of the veloc- 
ity fluctuations are actually independent of the details of the dissipative 
processes. The state depends only on the dimension of dissipative structures 
and the hierarchical relation index. 

5 Observable effects 

The nonlinearly evoloved IGM provides a coherent and uniform expanation of 
various observations, including the intermittence and non-Gaussianity of Lya 
transmitted fllux of quasar absorption espectrum [HJ US] ; Lya leaks of ab- 
sorption espectrum of high redshift objects[46j; the turbulence broadening [47J; 
Scaling relations between SZ effect and X-ray luminosity of X-ray clusters 
[381 HU S9] , baryon missing in collapsed halos [19]; etc. Here we address only 
two of them. 
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5.1 Non-Gaussianity of Log-Poisson hierarchy 



It has been shown that the SL scaling is yielded from the Log-Poisson hier- 
archy process, which is related to the so-called generalized scale covariance 
of the Navier-Stokes equations [50] . Therefore, the statistical behavior of 
the mass field of cosmic baryon matter should be given by the log-Poisson 
random multiplicative processes (RMP). 

The log-Poisson RMP assumes that, in the scale- free range, the variables 
5p r = p(x + r) — p(x), r = |r|, on different scales r are related from each 
other by a statistically hierarchy relation as 

\5pr 2 \ = W ri r 2 \5p ri \, (28) 

where 

W rrr2 = P m ( ri /r 2 )\ (29) 

which describes how the fluctuation \5p ri \ on the larger scale r 1 related to 
fluctuations \Sp r2 \ on the smaller scale r 2 . In eq.(29), m is a Poisson random 
variable with the PDF 



lm) 



exp(-\ rir2 )\™Jm\. (30) 



To insure the normalization (W rir2 ) = 1, where (...) is over m, the mean A rir2 
of the Poisson distribution should be 



\ rir2 = j[]n(n/r 2 ))/(l-/3). (31) 

It is enough to consider only \6p r \, as the distribution of positive and negative 
5p r is symmetric. 

The log-Poisson hierarchy given by eq.(29) depends only on the ratio 
r \/ r 2i which is obviously scale invariant. The hierarchy is determined by 
two dimensionless positive parameters: (3 and 7, describing, respectively, the 
intermittence and singularity of the random fields. Equation (29) relates Sp r 
on different scales by multiplying a random factor W, which generally yields 
a non-Gaussian field even if the field originally is Gaussian [5T] . 

The cascade from scale r\ to r 2 , and then to r 3 is identical to the cascade 
from r\ to r 3 . It is because W rirs = W rir2 W r2r3 = (3 N (ri/r 3 ) 7 , where N is 
again a Poisson random variable with A rir3 = A rir2 + A r2r3 = 7[ln(ri/r 3 )]/(l — 
0). Therefore, the log-Poisson hierarchy removes an arbitrariness in defining 
the steps of cascade from r± to r 2 or r 2 to r 3 . The log-Poisson hierarchy is 
discrete in terms of the discrete random number m. However, the scale r is 
infinitely divisible. Namely, there is no lower limit on the difference r% — r 2 . 
It can be infinitesimal, and the hierarchical process is of infinite divisibility. 
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With the log-Poisson hierarchy eqs.(28)-(31), the density field should have 
the non-Gaussianities, which can be described by the structure function of 
density field S p (r) = (\Sp r \ p ) as follows [52] : 

1. Intermittent exponent. The structure function should be a function of 
power law of r as 

S p (r) oc r^ p \ (32) 
and the intermittent exponent £(p) is given by 

£(p) = -7b -(l-/3 p )/(l (33) 

This is actually the SL scaling formula for density field. 

2. Hierarchical relation. Defining 

F p (r) = S p+1 (r)/S p (r), (34) 
the log-Poisson RMP gives hierarchical relation as 

(35) 



F r 



oo \ 



Fp+i(r) 
F r 



■oo \ 



K p = -7 ^ ^ T ^ > . (37) 



ri,ri — / x „P \ / x „P \ ' v 



Relation (35) is invariant with respect to a translation in p. 

3. Moments. The ratio between high order and 2nd order moments, 
(5p?P)/(5p 2 r ) p should satisfy 

(Sp 2p ) 

In = KAnr + const (36) 

with 

P (i - /3 2 ) - (i - p 2 n 

1-/3 

4. Scale-scale correlation. The so-called scale-scale correlation is defined 

as 

{5fi 1 ){5fi 2 )' 

Obviously, for a Gaussian field, C p f r2 = 1. If the ratio r 2 /r 1 is fixed, the 
log-Poisson model predicts the scale-scale correlation to be 

C™ n = Bin/nyf^ 2 ^, (39) 

where the coefficient B(r 2 /ri) depends only on the ratio Tijr\, as the log- 
Poisson model is invariant of the dilation. 

It should be emphasized again that all above-mentioned non-Gaussianities 
depend only on two dimensionless parameters /3 and 7. It has been revealed 
that the fields of Lya transmitted flux fluctuations of observed high resolution 
and high S/N Lya absorption spectra of quasars at moderate redshift z ~ 2 
- 3, and high redshift z ~ 5 - 6 do show all the log-Poisson non-Gaussianities 
[531 Ell- 
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5.2 Turbulence pressure and baryon missing 

As shown in §2.2, the Jeans length Aj = lnjkj is Aj = c s y/ir/Gp, c s being 
the sound speed. The effect of turbulent motions on the Jeans length can be 
estimated by replacing c s by an effective sound speed 

< eff = c 2 s + \(v 2 ) (40) 

where (v 2 ) is the rms velocity dispersion due to turbulent motion [551 156] . 
Namely, the random velocity field of turbulence contributes to an extra pres- 
sure, the turbulence pressure, p tu b = p(v 2 ) to resist the gravitational collaps- 
ing. 

If the collapsed halo is on scale R, the gravitational collapsing will not be 
affected by the rms velocity dispersion on scales larger than R, the fluctuation 
of velocity on the scales k < 2ir/Rdo not contribute to the turbulent pressure 
in terms of resist gravitational collapsing on scales that larger than R. Thus, 
the effect of turbulent pressure on gravitational clustering should be given 
by |5ZI 

/"^max 

Ptur(A*) = / E(k)dk, (41) 

JkR 

where E(k) is the power spectrum of kinetic energy of the turbulence, and 
k R = 2n/R. Wavenumber /c max = 27r//diss corresponds to the minimal scale 
/diss below which the turbulence decays due to energy dissipation or virializa- 
tion. The turbulent pressure is dynamical, not thermal. It can be as strong 
as a thermal pressure of temperature 10 5-6 K, while the temperature of the 
IGM is still remained in the range 10 4-5 K. 

Obviously, the bulk velocity will not contribute to the (v 2 ), as it can be 
given by the motion of falling into a gravitational well. Therefore, we should 
pick up the component of velocity field that is regardless of gravity. It is 
mainly given by the vorticity of velocity field. In other words, the turbulent 
vorticity field can prevent the collapsing of baryon matter into the halos of 
dark matter. It leads to the deviation of the mass density fields of baryon 
matter from that of dark matter. 

With the samples produced by cosmological hydrodynamical simulation, 
we find that the turbulence pressure will play an important role of preventing 
the IGM clustering. In the IGM regions with 10 times mean overdensity, the 
turbulent pressure is equivalent to the thermal pressure of the baryon gas 
with temperature of about 1 x 10 5 K. The turbulent pressure is dynamical 
and non-thermal. When the turbulence dissipated, the kinetic energy of 
turbulence becomes the thermal energy. It yields the entropy in halos [38] . 
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The baryon fraction of dark matter halos is found to be a strong function 
of the halo mass M 2 oo, the mass within the range of total mass density 
larger than 200 pt- The baryon fraction is significantly decreasing with the 
mass of dark halos. For halos with mass M200 = 3 x 10 10 M , the baryon 
fraction is only about 0.1 of the cosmic mean of the concordance ACDM 
universe, and it further decreases to 0.03 for M 20 o = 10 10 M Q halos. This 
result is well consistent with observed baryon content of galaxies and galaxy 
groups. Therefore, the turbulent velocity field of baryon fluid would be the 
major factor leading to the absence of baryon matter in halos with mass less 
than 10 115 M Q . With this evolution picture, no pre-heating or reheating 
mechanisms are needed for explaining baryon missing. 

6 Conclusion 

The nonlinear evolution of cosmic baryon fluid is a central problem of cos- 
mology. The concordance ACDM model provides an effective paradigm of 
formation and evolution of cosmic structures. The dynamical behavior of the 
CDM component has been successfully revealed with analytical study and N- 
body numerical simulation. On the other hand, the dynamics of baryon fluid 
is less understood, especially in the nonlinear regime, because the involve- 
ment of complicated physical processes, such as cooling, heating, ionization, 
shock, turbulence etc. On the other hand, many observations are directly 
related to the nonlinear evolution of baryon fluid, such as the absorption 
spectra of quasars, galaxies and GRB at high redshift, Sunyaev-Zeldovich 
effect, the polarization of cosmic microwave background (CMB), X-ray back- 
ground etc. One of the biggest challenges in galaxy formation is also on the 
physical origin of the non-linear relation between baryon gas and low-mass 
dark halos. 

Thanks for the development of reliable cosmological hydrodynamical sim- 
ulation in the last decade, significant progress has been made in understand- 
ing the non-linear dynamics of cosmic baryon fluid. The turbulence behav- 
ior of cosmic baryon fluid in scale free regime is the key of understanding 
nonlinear evolution of the mass density and velocity fields. An important 
consequence of the turbulence is to yield statistical decoupling of the mass 
and velocity fields of baryon fluid from the underlying dark matter field. This 
result is extremely valuable to understand the low baryon fraction in galaxy 
clusters, the non-Gaussian features of quasar Lya transmitted flux, the Lya 
leaking of absorption spectra of high redshift objects, the shock heating of 
pre-virialized structures, and the strongly nonlinear dependence of the star 
formation efficiency on halo mass etc. Considering these complexities, the 
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nonlinear dynamics of the IGM also is essential for understanding the epoch 
of reionization and related observations. 
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